Quantum Monte Carlo study of MnO solid 
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Electronic structure of the manganese oxide solid is studied by the quantum Monte Carlo (QMC) methods. 
The trial wavefunctions are built using orbitals from unrestricted Hartree-Fock and Density Functional Theory, 
and the electron-electron correlation is recovered by the fixed-node QMC. The correlation effects are significant 
and QMC estimations of the gap and cohesion show a very good agreement with experiment. Comparison 
with hybrid functional results points out the importance of the exact exchange for improvement of the Density 
Functional description of transition metal oxide systems. 

PACS numbers: 71.27.4*, 72.80.Ga, 71.20.-b, 71.15.Nc 



Transition metal compounds and transition metal ox- 
ides (TMOs) in particular belong to the most complex and 
important types of solid materials. TMOs exhibit a mul- 
titude of collective effects such as ferro-, ferri- and anti- 
ferromagnetism, ferroelectricity, superconductivity, etc 00]. 
The electronic structure of these systems poses a real chal- 
lenge both for theory and experiment and TMOs have re- 
mained on a forefront of condensed matter research for 
decades 0HEE0HI3II3EI1]. Among the TMO solids, 
MnO and NiO have become paradigmatic examples of strong 
electron-electron correlation effects and antiferromagnetic or- 
dering. Both MnO and NiO should nominally exhibit par- 
tially filled levels which would lead to a metallic ground state, 
however, experiments revealed that these systems are actually 
antiferromagnetic insulators with large gaps. The qualitative 
picture of electronic structure is explained either by a Mott- 
Hubbard mechanism implying that the gap results from a large 
local Coulomb repulsion in doubly occupied d- states, or by 
a charge transfer from transition metal d- states to oxygen p 
levels with the gap having p — > d like character or a by a 
combined regime in between 0]. The localized d subshells 
exhibit unpaired spins which at low temperatures order into an 
antiferromagnetic AF II insulator with alternating spin (111) 
planes in cubic rocksalt structure. 

The electronic structure of MnO and NiO have been stud- 
ied by a number of theoretical approaches, most notably by 
the spin-polarized Density Functional Theory (DFT) with sev- 
eral types of functionals. Augmented spherical wave lo- 
cal density approximation calculations |5|] and subsequent 
works |01 explained the stability of the AF II ordering in 
NiO and MnO. Very recently, the rhombohedral distortion 
of MnO was successfully described by accurate DFT cal- 
culations |4|]. However, the commonly used DFT function- 
als are less reliable for predicting other key properties of 
strongly correlated systems. For example, DFT band gaps 
can be underestimated by a factor of 3-10 or even absent, 
leading thus to false metallic states. Therefore a variety^of 
DFT modifications such as self-interaction correction 
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orbital-polarization corrections [9], on-site Coulomb interac- 
tion (LDA+U approach) |[l(J, [ill] and others, have been sug- 
gested to improve the description of gaps and other electronic 
properties. On the other hand, Towler et al lfl2h studied MnO 
and NiO with unrestricted Hartree-Fock theory (UHF). UHF 
is seemingly a poor method since it neglects the electron cor- 
relation completely. Nevertheless, for transition elements the 
exchange, which is exact in the HF theory, is at least as im- 
portant as correlation, especially for metallic ions with an ef- 
fective d— subshells occupation close to half-filling. The UHF 
results confirmed the crucial role of exchange in TMOs and 
provided a complementary picture to DFT with overestimated 
gaps and underestimated cohesion, but also with the correct 
AF order, magnetic moments within 10% from experiments 
and reasonably accurate lattice constants. Moreover, unlike 
DFT approaches which predict insulator only for the AF II 
ground state, UHF keeps the gap open also for the ferromag- 
netic or any spin-disordered phases. This agrees with exper- 
iment which shows that MnO is an insulator well above the 
Neel temperature T/v ~ 118K since the spin ordering Mn- 
O-Mn superexchange mechanism is very weak and hardly af- 
fected by spin flips on localized ions. 

In this letter, we take a fresh look on the MnO system in 
the framework of many-body quantum Monte Carlo method 
(for example, Ref. [13] and references therein). Our aim is 
to understand and quantify the impact of explicit treatment of 
both exact exchange and correlation on the key properties such 
as cohesion and band gap. In sp systems the QMC method 
was very successful in capturing electron correlation effects 
for a large number of valence electrons 111 311 . It is therefore 
both important and interesting to test the QMC performance 
on challenging problems such as TMOs, extending the previ- 
ous attempts to apply QMC to NiO system II 1 4 1 1 511 . 

We employ variational Monte Carlo and the diffusion 
Monte Carlo (DMC) methods. The trial/variational wavef unc- 
tion is expressed as a product of Slater determinants of single- 
particle spin-up and spin-down orbitals ({^a},{^/3}) multi- 
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plied by a Jastrow correlation factor I16LI17I1 . 
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where / corresponds to the ions, z, j to the electrons and 
rjj, T^ur^i to the distances. Similarly to the previous work 
B13[ Il4 Il8l1 the correlation function consists of a linear com- 
bination of electron-electron and electron-ion terms and in- 
volves 9 variational parameters. The DMC method is used 
to remove the major part of variational bias which is inher- 
ent to the variational methods. DMC is based on the property 
that the projection linv^oo exp(— tH)^t is proportional to 
the ground state for an arbitrary with the same symmetry 
and non-zero overlap fl9ll . The fermion sign problem caused 
by the antisymmetry of electrons is avoided by the commonly 
used fixed-node approximation. 

In our calculations the core electrons are eliminated by 
pseudopotentials. For transition elements the choice of accu- 
rate pseudopotentials is far from trivial. In many TMO pseu- 
dopotential calculations with plane wave basis sets only the 
outermost d and s states are included into the valence space. 
However, earlier work [ 21 ] indicated that the errors from elim- 
ination of the so-called semicore states, such as 3s and 3p in 
the 3d transition series, could be significant. To elucidate this 
aspect, Table I. provides a comparison of the s — > d excitation 
energies for the Mn atom as obtained by the DMC method us- 
ing scalar relativistic large-core (Ar) 1 22] and small-core (Ne) 
pseudopotentials 1E3I1 . It is evident that inclusion of 3s and 
3p states into the valence space is important. The large-core 
(Ar) pseudopotentials show errors of the order w 0.5-0.6 eV 
while the Ne-core pseudopotentials are much more accurate 
with typical errors ^0.1 eV, what is comparable to the bias 
from the fixed-node approximation 111 311 . The physical reason 
of such behavior is well-known and stems from the spatial dis- 
tribution of the 3s, 3p electrons which is similar to the one of 
3d electrons since they occupy the same principal shell. The 
price which one has to pay for using small core is, of course, 
significant. The magnitude of total energy, which is one of 
the key measures of QMC computational demands, increases 
roughly eightfold. It is also useful to notice that the DFT re- 
sults for these excitations differ from experiments by 0.6-1 eV, 
clearly indicating the large contributions from exchange and 
correlation (Tab. I.). 

For the MnO solid we first carried out calculations with the 
spin-unrestricted Hartree-Fock and DFT (B3LYP and PW86) 
methods and Ne-core pseudopotentials using the CRYS- 
TAL98/03 packages B25I1 . The orbitals were expanded in 
gaussian basis sets with (12s, 12p, 7d) gaussians contracted to 
[3s, 3p, 2d] and (8s, 8p, Id) contracted to [4s, 4p, Id] for Mn 
and O atoms, respectively. Figure 1 shows the band structure 
of MnO solid which is obtained from UHF (a), B3LYP (b), 
PW86 (c) methods. Note that B3LYP hybrid functional con- 
tains 20% of the Hartree-Fock exchange so that it "interpo- 
lates" between the exact HF exchange and the effective local 



TABLE I: The excitation energies [eV] s — > d of the Mn atom 
as calculated by UHF, DFT (BLYP and B3LYP) and DMC meth- 
ods. We used the Ne-core scalar relativistic pseudopotentials except 
for the DMC(Ar) calculations which employed the Ar-core pseu- 
dopotentials. In the last column are the experimental energies. Our 
DFT/BLYP results are close to similar calculations done previously, 
seeRef. IH. 

UHF BLYP B3LYP DMC(Ar) DMC(Ne) Exp. 



3.5 
9.1 



1.2 
4.4 



1.6 
5.2 



1.6(1) 
5.2(1) 



2.2(1) 
5.8(1) 



2.1 
5.6 



DFT exchange limits 12611 and it often provides an improved 
picture of excitations both in molecules and solids. 

In QMC the MnO solid is represented by a supercell with 
periodic boundary conditions. This way of simulating the ac- 
tual solid involves finite size errors which scale as 1/N where 
N is the number of atoms in the supercell |23,l231- The fi- 
nite size errors affect mainly the estimation of cohesive energy 
where one needs to calculate the energy per primitive cell vs. 
isolated atoms. In order to filter out the finite size bias we have 
carried out VMC calculations of supercells with 8, 12, 16, 20 
and 24 atom/supercell. The most accurate and extensive fixed- 
node DMC calculations were carried out with B3LYP orbitals 
for 16 and 20 atoms in the supercell. The cohesive energy 
obtained by the DMC method shows an excellent agreement 
with experiment (Tab. II). To evaluate the impact of the cor- 
relation on the gap we have estimated the energy of the T — > 
B excitation by an exciton calculation B18I1 . The QMC result 
for excitation energy is less perfect with the difference from 
experiment being a fraction of an eV. This clearly shows that 
the wavefunction and corresponding fixed-node error is larger 
for an excited state. Nevertheless, the differential energy gain 
for excited vs. ground state from correlation of « 8 eV is sub- 
stantial and demonstrates the importance of this effect both in 
qualitative and quantitative sense. Due to large computational 
demands, similar but statistically less precise DMC calcula- 
tions were carried out also with the UHF orbitals. While for 
the ground state the difference between the two sets of orbitals 
was marginal, the excited state with UHF orbitals appeared 
higher in energy approximately by w 1.5(0.5) eV indicating 
thus, not surprisingly, even larger fixed-node bias for the ex- 
cited state in the UHF approach. 



TABLE II: The MnO solid cohesive energy and B — > V excitation 
energy calculated by UHF, DFT and DMC methods compared with 
experiment. The determinantal part of the DMC wavefunction used 
the B3LYP one-particle orbitals. 



UHF PW86 B3LYP DMC Exp. 



E co h 6.03 
B -> T 13.5 



11.00 
1.2 



9.21 
4.0 



9.40(5) 9.50 
4.8(2) ^4.1 



It is quite encouraging that the fixed-node DMC with the 
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FIG. 1: The band structures of MnO obtained by (a) the unre- 
stricted Hartree-Fock method, (b) DFT/B3LYP functional, and (c) 
DFT/PW86 functional. The PW91 and PBE functionals provide es- 
sentially the same picture as the PW86 functional. The calculated 
excitation in QMC is indicated by an arrow on the UHF plot and 
the corresponding one-particle states are denoted by open and filled 
circles. 



simplest possible single-determinant wavefunction leads to a 
consistent and parameter- free description of the basic proper- 
ties of this strongly correlated system. An obvious question 
to ask is whether one-determinant is sufficiently accurate for 
an antiferromagnet since the wavefunction with different spin- 
up and spin-down orbitals is manifestly not an eigenfunction 
of the square of the total spin operator. In order to eliminate 
the spin contamination one would need to explore wavefunc- 
tion forms beyond the single-determinant S later- Jastrow, for 
example, generalized valence bond wavef unctions. However, 
since the actual mechanism is the Mn-O-Mn superexchange, 
one can expect the resulting effect to be small and most prob- 
ably undetectable within our error bars (our calculations with 
wavefunctions beyond the single UHF determinant were not 
conclusive.) 

It is interesting to revisit now the one-particle results and 
provide some feedback from our QMC calculations. The anal- 
ysis of orbitals indicates that the nature of the top valence 
bands is rather similar in all approaches with both p and d 
states having significant weights in these states across the Bril- 
louin zone (see Fig. 2). This is clear also from the Mullikan 
population analysis which shows effective magnetic moments 
on Mn atoms in UHF, B3LYP and PW86/PW9 1/PBE meth- 
ods to be 4.92, 4.84 and 4.78/iB, respectively; these values 
are quite close to each other and border the range of experi- 
mental estimates of 4. 58 -4.78 ps- 

The bottom of the conduction band is free-electron-like 
T state with significant amplitudes from atomic 0(3 s) and 
Mn(4s) orbitals and it is this state which is responsible for the 
DFT gap closing in ferromagnetic or spin-disordered phases. 
For the ferromagnetic phase B3LYP exhibits a gap of w 2.4 
eV and it is straightforward to check that by decreasing the 
weight of exact exchange the gap decreases. For example, 



UHF DFT/B3LYP 




FIG. 2: Isosurfaces of UHF and DFT/B3LYP (blue/dark region is 
positive while yellow/light is negative). The valence B states are 
plotted for four atoms while the V states are plotted for a supercell so 
that their conducting character is visible. The valence states have sig- 
nificant weights both from Mn d states and O p states. Note small in- 
crease/decrease of d/p orbital amplitudes from UHF to DFT/B3LYP. 
The conduction state is much more delocalized and is composed 
from Mn 4s and O 3s atomic states. 

with 10% of exact exchange in B3LYP the gap lowers to « 
1.2 eV. The functionals without the exact exchange, such as 
PW86/PW9 1/PBE, lead to ferromagnetic metals due to the 
overlaps with the uppermost valence bands and subsequent 
rehybridization of states around the Fermi level. This is the 
"usual" DFT band gap problem which favors smooth and ex- 
tended states at the expense of the more localized ones which 
are stabilized by the exact (nonlocal) exchange and correlation 
effects. 

The results presented here therefore provide quite sugges- 
tive insights into the problem of a simple one-particle model 
appropriate for TMOs. In the Mott-Hubbard picture the ori- 
gin of the band gap is the large on-site repulsion between 
the d electrons which basically relates the gap to the two- site 
d n d n — > d n ~ 1 d n + 1 type of excitation. On the other hand, the 
charge transfer favors the ionic picture of Mn ++ __ with the 
oxygen p states at the top of the valence band the gap given by 
the p — > d excitation energy. Our QMC results quantify that 
UHF produces qualitatively correct wavefunctions although 
biased towards the charge transfer limit, especially for excited 
states. The hybrid B3LYP functional provides more balanced 
zero-order theory, arguably better and less biased than DFT 
non-hybrid functionals. This view is supported also by the 
orbital analysis proposed by Brandow dtj, and, in effect, 
also by introduction of on-site terms which restore some of 
the Hartree-Fock character. To a certain extent, the hybrid 
functional alleviates the DFT biases by eliminating part of the 
self-interaction and by introducing the exchange "pull-down" 
attraction resulting in lower energies of localized states. Us- 
ing such orbitals in QMC correlated framework enables us to 
obtain results which are close to experiment without any ad- 
ditional parameters. 
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It is tempting to consider the usefulness of a hybrid func- 
tional which would be, however, derived from a fundamental 
theory instead of a fit to a testing set of molecules underly- 
ing B3LYP B26I1 . Although such approach would not fix all 
the deficiencies of the approximate DFT functionals it could 
serve, for example, as a cost effective method for providing 
more appropriate sets of one-particle orbitals for building ac- 
curate wavef unctions. In fact, for molecular systems such as 
TiO and MnO we were able to directly optimize the weight 
of the exact exchange within a QMC framework. By vary- 
ing the weight of the exact exchange in the functional and by 
iterating fixed-node DMC calculations we found the best set 
of one-particle orbitals which provided the lowest fixed-node 
energy in DMC 12711 : for solids such calculations might be 
possible in the near future. 

In conclusion, we have carried out calculations of the MnO 
solid in the variational and diffusion Monte Carlo methods. 
We have evaluated the cohesive and excitation energies which 
show excellent agreement with experiment. We have pointed 
out the necessity of using high accuracy Ne-core pseudopo- 
tentials for Mn because of semicore 3s and 3p states. The 
results clearly show a crucial role of both exchange and cor- 
relation and their accurate description not only for MnO but 
obviously for other transition metal oxide systems as well. 
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